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ABSTRACT 


i 

An  analytical  model  is  presented  for  estimating  the  strearowise  length  of 
the  steady-state  separation  bubble  induced  by  the  interaction  of  a  shock 
moving  with  constant  speed  over  a  thin  thermal  (high  sound  speed)  layer.  In 
the  case  of  a  thermal  layer  of  semi-infinite  extent,  the  steady  solution 
represents  the  asymptotic  limit  of  an  initially  unsteady  interaction.  The 
streamwise  length  of  the  separation  bubble  is  found  by  equating  mass  flow  into 
the  bubble  to  the  mass  flow  removed  by  the  wall  boundary  layer.  Numerical 
results  are  presented  for  a  variety  of  shock  speeds,  thermal  layer  heights, 
and  thermal  layer  sound  speeds.  The  present  study  indicates  the  importance  of 
viscous  effects  on  precursor  development  and  provides  a  method  for  modifying 
the  wall  boundary  condition  in  inviscid  numerical  codes  in  order  to  include 
the  effect  of  the  wall  boundary  layer  on  the  inviscid  flow  field. 
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I.  INTRODUCTION 


The  interaction  between  a  moving  shock  wave  and  a  thin  thermal  (high 
sound  speed)  layer  results  in  a  precursor  wave  system,  as  first  discussed  by 
Hess.^  Experimental  studies  of  this  interaction  are  reported  in  Refs.  2-^J. 
Numerical  code  calculations,  which  assume  inviscid  flow,  are  presented  in 
Refs.  5-7. 

In  Ref.  8,  a  simple  analytical  model  was  presented  that  neglected  viscous 
effects  and  that  was  in  approximate  agreement  with  the  inviscid  code  calcula¬ 
tions  in  Ref.  5.  The  inviscid  model  and  the  inviscid  code  calculations  indi¬ 
cate  that  the  scale  of  the  interaction  increases  linearly  with  time.  In  Ref. 
8,  it  was  pointed  out  that  viscous  effects  are  expected  to  cause  the  scale  to 
approach  an  asymptotic  limit.  In  this  limit,  the  flow  is  steady  in  shock 
fixed  coordinates.  The  asymptotic  scale  is  estimated  in  the  present  study  by 
consideration  of  wall  boundary  layer  effects.  The  time  to  reach  steady  state 
is  also  estimated.  Finally,  the  application  of  the  present  results  to  invis¬ 
cid  codes  is  discussed.  Symbols  are  defined  in  the  Appendix. 
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1 1 .  THEORY 


Consider  a  normal  shock  wave  that  moves  with  uniform  velocity  u  over  a 

3 

thin  thermal  layer  of  semi-infinite  extent.  Flow  conditions  prior  to  the 
start  of  the  interaction  are  indicated  in  Fig.  1,  using  a  laboratory  station¬ 
ary  coordinate  system.  In  this  coordinate  system,  flow  velocities  are  denoted 
by  a  superscript  bar.  The  thermal  layer  is  characterized  by  a  speed  of  sound 
a^,  which  is  higher  than  the  adjacent  fluid  value  a-j. 

For  cases  other  than  a^/a^  -  1  small  and/or  ug/a^  -  1  small,  a  shock  is 
generated  in  the  thermal  layer  that  moves  faster  than  the  incident  shock  and 
thereby  creates  a  "precursor"  flow  field.1’®  The  precursor  flow  field  at 
small  times,  when  viscous  effects  are  negligible,  is  illustrated  in  Fig.  2, 
using  a  coordinate  system  wherein  the  incident  shock  is  stationary.  Veloci¬ 
ties  in  this  coordinate  system  are  denoted  by  the  unbarred  symbol  u.  The 
total  head  of  the  thermal  layer  fluid  in  region  A  is  less  than  the  static 
pressure  in  region  2.  Consequently,  the  high  pressure  gas  from  region  2 
expands  into  region  3  and  drives  a  shock  in  the  thermal  layer.  The  leading 
edge  of  the  driver  gas  is  denoted  Xj,  and  the  thermal  layer  shock  location  is 
denoted  xst.  The  quantities  xgp  and  xg  denote  a  rear  stagnation  point  and 
incident  shock  location,  respectively.  The  thermal  layer  gas  collects  in  the 

form  of  a  bubble  beneath  the  incident  shock.  Inviscid  numerical  code  calcula¬ 
te 

tions-'  indicate  that  the  interface  velocity  u^  and  the  thermal  layer  shock 
velocity  usfc  are  equal  and  remain  constant  with  time.  As  a  consequence,  the 
separation  distance  xgt  -  x^  remains  constant.  An  analytic  model  of  the 
precursor  flow  field  was  presented  in  Ref.  8  by  analogy  with  a  low  pressure 
shock  tube  flow  wherein  the  separation  distance  between  shock  and  contact 
surface  remains  constant. 

Inviscid  solutions  for  the  precursor  flow  field  indicate  that  the  scale 
of  the  interaction  grows  linearly  with  time.  As  noted  in  Ref.  8,  viscous 
effects  will  cause  the  scale  to  approach  an  asymptotic  limit.  In  this  limit, 
the  rate  of  mass  flow  into  the  interaction  region  equals  the  rate  of  mass  flow 
out  of  the  interaction  region  because  of  the  wall  boundary  layer.  This  steady 
state  asymptotic  limit  is  illustrated  in  Fig.  3t  using  an  incident  shock  fixed 
coordinate  system.  In  this  coordinate  system,  the  wall  moves  with  velocity 
Uy  =  u1  =  u  .  The  wall  boundary  layer  may  be  viewed  as  aspirating  mass  from 
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Fig.  1.  Initial  Conditions  in  Laboratory  Stationary  Coordinate  System 


'////////. 


Fig.  2.  Flow  Conditions  in  Incident  Shock  Stationary  Coordinate  System  at 
Early  Times.  Viscous  effects  are  negligible,  and  the  scale  of  the 
interaction  grows  linearly  with  time. 
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Fig.  3.  Flow  Conditions  in  Incident  Shock  Stationary  Coordinate  System  at 
Late  Times.  Wall  boundary  layer  flow  is  illustrated.  The  scale  of 
the  interaction  is  fixed,  and  the  flow  is  steady. 
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the  interaction  region.  The  limiting  streamline  in  Fig.  3  separates  the  free 
stream  mass  flow,  which  enters  the  bubble  region,  from  the  free  stream  mass 
flow,  which  does  not  enter  the  bubble  region. 

The  asymptotic  length  of  the  interaction  region,  denoted  L,  is  estimated 
in  the  following  sections  by  equating  the  mass  flow  into  and  out  of  the  bubble 
region . 

The  rate  at  which  mass  enters  the  bubble  region  (per  unit  width)  may  be 
written 

“in  =  p4u4h4  +  p1u1h1  { 

—  PjjtijjhjjB 


where,  noting  u1  =  u^ 


B  =  1  +  [  (P1h1  ] 

Here,  hjj  is  the  thickness  of  the  thermal  layer,  and  h^  is  the  thickness  of  the 
free  stream  layer  that  enters  the  bubble  (Fig.  3).  The  quantity  B  is  the 
ratio  of  total  inflow  to  the  inflow  associated  with  the  thermal  layer.  The 
quantity  Bh^  may  be  viewed  as  an  effective  height  of  the  thermal  layer. 

The  wall  boundary  layer  in  Fig.  3  is  reproduced  In  Fig.  4a.  The  solution 
of  this  boundary  layer  flow  is  complicated  by  the  presence  of  the  contact 
surface  (which  separates  free  stream  particles  from  thermal  layer  particles) 
and  by  nonuniform  boundary  layer  edge  conditions.  For  present  purposes,  the 
boundary  layer  effect  will  be  estimated  by  considering  the  limit  Ug/u^,  -»  o  for 
a  boundary  layer  with  uniform  edge  conditions  (Fig.  4b).  This  approach  is 
accurate  for  the  case  of  strong  incident  waves  (wherein  ug/uw  and  u^/uw  are 
near  zero)  and  provides  an  average  edge  condition  for  the  case  of  weak  inci¬ 
dent  waves.  Distance  from  the  shock  wave  in  the  thermal  layer  is  denoted  t. 

The  excess  mass  flow  in  the  boundary  layer  at  any  streamwise  station  is 

6 

“bl  =  /  (flU  _  peue)dy  (2) 

o 

where  6  is  the  local  boundary  layer  thickness.  The  excess  mass  flow  in  the 

boundary  layer  may  be  expressed  in  the  form 
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in  Incident  Shock  ^ixed  boundary 

;c"cy  r^resen^ation  actual ^  ^  ^ry  layer. 

(b)  Assumption  of  uniior 
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* 

where  6  is  the  boundary  layer  displacement  thickness 


* 

5 


)dy 


(3b) 


The  vertical  velocity  at  the  outer  edge  of  the  boundary  layer  is,  from  conti¬ 
nuity  considerations, 


p  v  = 
e  e 


dmbl  _  (  dfi*\ 

d$  peue^~  d£  ' 


(4) 


*  Q 

Both  6  and  v„  are  negative  for  shock  induced  boundary  layers.  Solutions 

*  c 

for  6  and  vg  are  presented  in  Ref.  9  for  the  case  of  uniform  flow  external  to 
the  boundary  layer.  If  it  is  assumed  that  the  density  and  viscosity  within 

the  boundary  layer  are  characterized  by  the  wall  values,  p  and  p  ,  the  inte- 

ft  w  w 

gral  solution  for  6  and  vg  in  Ref.  9  indicates  that  for  t 


-^1  -2  —  =  I.™  (V/2 

p  u  £  p  u  vu  Z' 

w  w  w  w  w' 


(5a) 


and,  for  t  t 


m, 


bl 


p  u  U  -  f.  )' 
w  w  ^o 


1/5 


(5b) 


where  Eqs.  (5a)  and  (5b)  apply  for  laminar  and  turbulent  boundary  layers, 
respectively.  Here,  refers  to  the  effective  origin  of  the  turbulent  bound¬ 
ary  layer,  and  refers  to  the  location  of  the  transition  from  a  laminar  to  a 
turbulent  boundary  layer  (Fig.  5).  The  transition  is  assumed  to  be  instan¬ 
taneous  and  is  assumed  to  occur  at  a  known  Reynolds  number  defined  by 


Vt 


=  Re. 


(6) 


If  it  is  assumed  that  m&1  is  continuous  at  the  quantities  and  are 
related  by 
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Fig.  5 


Boundary  Layer  Development  with  Transition  at  and  with  Effective 
Origin  of  Turbulent  Boundary  Layer  at 
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(7) 


Co  =  *t<1  -  20.28  Ret_3/8) 

Equations  (5)- (7)  define  the  excess  mass  flow  in  the  boundary  layer.  Equating 
and  m^n  provides  the  length  of  the  separation  bubble.  The  result  is, 
noting  uw  =  Ui,  =  u1f 


UlBh4  Bh4 


=  ( 


V‘,M) 

ill  j  v 

w 


1.404- 


LL  <  «t 


(8a) 


1/4 


4  "  4  ,v 4 


Bh, 


,  5/4  1/4 

°n/pu  vi 

=  ^0.1264^  W 


LL  >  «t 


(8b) 


where  and  Lj.  denote  the  length  of  the  laminar  and  the  laminar  plus  turbu¬ 
lent  bubble,  respectively.  The  wall  temperature  Tw  may  be  assumed  to  remain 
at  its  initial  value  T^.^  We  also  assume  that  regions  1  and  4  contain  the 
same  type  of  gas  (e.g.,  N2/N2  rather  than  N2/He) ,  and  introduce  the  expres¬ 
sions 


fi  fi  L.  Ei 

°u  '  Tit  '  »e 

q.W 

V  -  T 


(9a) 

(9b) 


Ms  --  U,/«, 


H  2  10"5  «  Bh.  a ,/v, 


into  Eq.  (8).  For  laminar  flow, 


J.  4.  .  0.5073  x  105 
H  Bh4=  pg/p1 


Rs  (VV 


1+UI 


1_  4  _  0.5073  x  1010  *s 

H2  4  Ret  (T^/T,)20"-10* 


(9c) 

(9d) 


(10a) 


(10b) 
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For  turbulent  flow, 


1  VJo  236.0  ,1/lt  ,Tu(1^)/4 

h1/4  Bh,  '  pe/P)  '"s'  1?^ J 


(11a) 


HJo  10'5  *  Wl^  f.  20.28, 

%  =  (P./P,».  0  '  «et3'8) 


(11b) 


Finally,  we  assume  that  the  bubble  pressure  pe  is  equal  to  the  stagnation 
pressure  in  region  5  of  Fig.  3.  For  an  ideal  gas, 


£e  f,  T*  ~  1  .^W11 

P,  P„  =  1  *  2  V 


2  4- 


(T„  ♦  DMj  V(V1}  /  Tj|  ♦  1  1/(Yi|-1) 

=  2  5  w2  J 

2y4MJ1  _  Yjj  ♦  1 


MH  1  1 


m4  >  1 


(12a) 


(12b) 


Equations  (10)-(12)  can  be  evaluated  by  specification  of  Ms,  Ti+/T1 1  H,  Ret, 
Yjj,  and  u.  Typically, 


and  for  air 


Ret  =  105 


=  7/5 


(14a) 


u  =  0.7 


(14b) 


-1  =  2.274  *  105  p  (222)  ’  cm'1 

i  M 


(14c) 


where  p1  and  Tj  are  in  atmospheres  and  K,  respectively.  The  quantity  H  is  in 
the  form  of  a  Reynolds  number  and  is  of  order  1  when  Bhjj  is  of  the  order  of  a 
centimeter. 


Numerical  results  are  given  in  Table  1  and  Fig.  6  for  1 .5  <  M  <  6.0,  1.0 
<  ajj/ai  £  3.4,  and  0.1  <  H  <  10.  For  each  Ms,  the  smallest  value  of  ajj/a-j  in 
Table  1  corresponds  to  the  onset  of  the  separation  bubble.  When  H  <<  1,  the 
boundary  layer  is  wholly  laminar  and  H  L^/CBh^)  is  of  order  10J  to  10  . 

When  H  =  0.1,  the  boundary  layer  is  laminar  for  the  larger  values  of  a^/ai  and 
is  laminar/turbulent  for  the  smaller  values  of  a^/a^.  In  the  latter  case, 
Lq./(Bhjj)  is  of  order  101  to  10^,  with  the  smaller  values  of  L^/ 

(Bh|j)  occurring  at  the  larger  value  of  Ms.  When  Bh^  is  of  the  order  of  1  cm,  H 
is  of  order  1  and  the  boundary  layer  is  laminar/turbulent .  The  parameter 
L-p/(Bh^)  may  be  interpreted  as  defining  the  precursor  length  Lp  in 
centimeters.  The  latter  case  then  has  values  in  the  range  8.5  <  Lp,  cm  <  412 
for  Ms  <  6.0  in  Table  2.  With  increase  of  H  beyond  10,  the  boundary  layer  is 
predominately  turbulent  and  the  parameter  H_1/^  Lp/(Bh^)  becomes  independent 
of  H. 


13 


Table  1.  Asymptotic  Length  of  Thermal  Precursor  for  Case  y 
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Table  1.  Asymptotic  Length  of  Thermal  Precursor  for  Case  =  7/5,  u»  =  0. 7  and  Re 
Laminar  length  LL  applies  when  0.000  Indicated  for  Lj.  (Cont'd) 
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Table  2.  Characteristic  Time  to  Reach  Steady  State 
for  Case  y1  =  =  7/5 


M 

s 

a4 

a1 

u. 

i 

a4 

Ref.  8 

u  t 
s  c 

ll,t 

Eq.  18 

1.500 

1 .232 

1.461 

5.001 

1.500 

1  .-400 

1.405 

3.21 2 

1.500 

1 .800 

1.252 

1.990 

1.500 

2.200 

1.110 

1.592 

1.500 

2.600 

0.989 

1 .400 

1.500 

3.000 

0.886 

1.295 

1.500 

3.400 

0.800 

1.230 

2.000 

1.135 

1.975 

8.277 

2.000 

1.200 

1.950 

5.882 

2.000 

1.400 

1.851 

3-382 

2.000 

1.800 

1.632 

2.133 

2.000 

2.200 

1.438 

1.719 

2.000 

2.600 

1.278 

1.512 

2.000 

3.000 

1.147 

1.388 

2.000 

3.400 

1.037 

1 .311 

2.500 

1.101 

2.478 

10.952 

2.500 

1 .200 

2.429 

6.027 

2.500 

1 .400 

2.294 

3.513 

2.500 

1.800 

2.006 

2.251 

2.500 

2.200 

1.760 

1.822 

2.500 

2.600 

1.560 

1.607 

2.500 

3.000 

1.398 

1 .476 

2.500 

3.400 

1.264 

1.391 
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III.  TIME  TO  REACH  STEADY  STATE 


We  now  estimate  the  time  to  reach  steady  state  for  the  interaction  of  a 
shock  with  a  thin  semi-infinite  thermal  layer  (Fig.  II.  The  bubble  growth  rate 
is  most  rapid  at  small  times,  when  the  interaction  is  essentially  inviscid. 
It  is  this  small  time  growth  rate  that  is  used  to  characterize  the  overall 
rate  of  bubble  growth. 

The  inviscid  solution  of  Ref.  8  provides  numerical  estimates  for  the 
unsteady  streamwise  separation  between  the  interface  and  the  incident  shock, 
namely , 

x .  -  x  =  t  (u.  -  u  )  (15) 

i  s  i  s 

where  u./a,,  is  tabulated  as  a  function  of  M  and  an /a,  in  Ref.  8.  The  time 
l  4  s  ^  1 

for  this  separation  distance  to  reach  a  value  LL  ^  is  denoted  tc  and  is  given 

by 

u  t  ,  a.  u,  . 

■f-  =  (Jrr  - ')  (,6) 

LL,T  M  a1  a4 
s 

when  Ll  ^  denotes  either  LL  or  Lp.  The  time  tc  is  assumed  herein  to  charac¬ 
terize  the  time  required  to  reach  the  asymptotic  steady  state  flow.  Note  that 
the  product  ustc  characterizes  the  distance  the  shock  moves  during  the  un¬ 
steady  portion  of  the  interaction. 

Numerical  results  for  ustc/L^  are  given  in  Table  2.  This  parameter  is 
of  order  10  for  values  of  a^/a^  that  correspond  to  the  onset  of  separation, 
and  approaches  a  value  of  order  1  with  increase  in  a^/a^. 
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IV.  APPLICATION  TO  NUMERICAL  CODES 


The  present  theory  suggests  a  method  by  which  wall  viscous  effects  can  be 
incorporated  into  inviscid  codes.  For  the  case  of  a  flow  that  is  steady  in 
shock  fixed  coordinates,  the  inviscid  code  wall  boundary  condition  vg  =  0  can 
be  replaced  by  the  condition,  from  Eqs.  (5a)  and  (5b), 


p  v 
e  e 


p1us 


al 


V  1/2 


5  <  *t 


hu  1/5 

'rl&r) 


t 


(17a) 


(17b) 


where 


«„  *  «t  -  2°-28  »n'3/8  ) 
and,  for  the  same  gas  in  regions  1  and  4, 


v,  1/2 

-M 


4-w 


V,  1/5 

-V) 


( 17c) 

(17d) 

( 17e) 

(17  f) 

(17g) 


The  quantity  Pe/Pi  is  evaluated  from  Eq.  (12a).  For  >  1,  the  origin  of 
the  c  coordinate  system  is  taken  to  be  at  the  shock  in  the  thermal  layer.  For 
<  1 ,  the  origin  is  taken  at  the  streamwise  location  where  a  prescribed 
pressure  increment  has  occurred. 
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Numerical  code  calculations  (either  viscous  codes  or  inviscid  codes  with 
prescribed  Pgve  *  0)  provide  a  basis  for  estimating  the  parameter  B,  which  was 
introduced  in  Eq.  (1).  Thus,  if  and  Lp  denote  steady-state  bubble  length 
values  obtained  from  numerical  codes,  the  corresponding  values  of  B  are 

Th  L  1/2 

8  =  2  LLl*t  (l8a) 

c  T,  L  -  i  4/5 

:  I  fj  ^~h "  ^  4  >  *t  (l8b) 

Hopefully,  B  is  a  weak  function  of  initial  conditions.  The  numerical  code 
calculations  can  also  be  used  to  obtain  improved  estimates  for  t£. 
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V.  CONCLUDING  REMARKS 


The  present  study  represents  a  first  attempt  to  analytically  estimate  the 
asymptotic  streamwise  length  of  the  separation  bubble  associated  with  the 
interaction  between  a  moving  shock  and  a  thin  stationary  thermal  layer. 
Parameters  and  qualitative  trends  are  identified.  As  previously  noted,  fur¬ 
ther  study  incorporating  numerical  code  computation  is  needed  to  provide  more 
accurate  estimates  of  the  asymptotic  bubble  length  and  of  the  time  to  reach 
steady  state  flow.  It  has  been  assumed  in  Eqs.  (9a)-(9d)  that  regions  1  and  4 
contain  the  same  gas  but  at  different  temperatures.  The  role  of  dissimilar 
gases  (e.g.,  N2/He)  also  requires  further  study. 

The  present  study  indicates  the  importance  of  viscous  effects  on  precur¬ 
sor  development  and  provides  a  method  for  modifying  the  wall  boundary  condi¬ 
tion  in  inviscid  numerical  codes  in  order  to  include  the  effect  of  the  wall 
boundary  layer  on  the  inviscid  flow.  Equations  ( 1 7a )  —  ( 1 7g >  are  approximately 
correct  for  cases  where  the  flow  is  steady  in  incident  shock  fixed  coordi¬ 
nates.  The  expressions  for  pgve  in  Eqs.  (17a)  and  (17b)  require  modification 
for  cases  where  the  flow  is  unsteady  in  shock  fixed  coordinates. 
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APPENDIX 


SYMBOLS 


a 

B 


H 

h1 

hH 

ll,lt 


in 

bl 


P 

t 

t 

u 

u 


°w 


lst» 

‘i* 


sp 


T 

6 

6 


V 

v 

C 


speed  of  sound 

ratio  of  total  flow  to  thermal  layer  flow  into  separation  bubble,  Eq. 
(1) 

normalized  thermal  layer  height,  10"^  Bh^  a-j/v^ 

thickness  of  free  stream  layer  that  enters  separation  bubble,  Fig.  3 
height  of  thermal  layer,  Fig.  3 

length  of  separation  bubble  for  laminar  and  laminar /turbulent  boundary 
layers,  respectively 

shock  Mach  number  in  wall  stationary  coordinate  system,  us/a1 

Mach  number  of  flow  in  thermal  layer  in  incident  shock  stationary 
coordinates,  ujj/a^  =  (a i /a^ )  Mg 

mass  rate  of  flow  into  separation  bubble,  Eq.  (1) 

excess  mass  flow  rate  in  wall  boundary  layer,  Eq.  (2) 

static  pressure 

time 

characteristic  time  to  establish  steady  flow,  Eq.  (16) 

streamwise  velocity  in  wall  stationary  coordinate  system,  Fig.  1 

streamwise  velocity  in  incident  shock  stationary  coordinate  system, 
Fig.  2. 

wall  velocity  in  incident  shock  stationary  coordinate  systems.  Fig.  3 
velocity  in  +y  direction  at  edge  of  boundary  layer 

streamwise  distance  measured  from  leading  edge  of  thermal  layer,  Fig.  1 

location  of  shock  in  thermal  layer,  interface,  incident  shock,  and 
stagnation  point,  respectively,  Fig.  2. 

distance  normal  to  wall 
ratio  of  specific  heats  Cp/cv 
boundary  layer  thickness 

boundary  layer  displacement  thickness,  Eq.  (3b) 
viscosity 

kinematic  viscosity,  y/p 

distance  measured  from  shock  in  thermal  layer,  xst-x 
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APPENDIX  (continued) 


p  density 

u  viscosity  exponent,  m  -  Tu 

Subscripts 

1,2, 

3,M,  flow  regions,  Fig.  2 

5 

e  edge  of  boundary  layer 

w  property  at  wall 

Superscripts 

M,u  barred  quantities  are  in  wall  fixed  coordinate  system 

M,u  unbarred  quantities  are  in  incident  shock  fixed  coordinate  system 


